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Veronicastrum sibiricum var. zuccarinii (Plantaginaceae) is an endangered perennial herb native to 
semi-natural grasslands in southern Korea and Kyushu, Japan. To elucidate the genetic status of V sibiri¬ 
cum var. zuccarinii in the Aso region of Kyushu, in which relatively large populations still remain, we 
analyzed seven populations using five novel microsatellite markers. Each population of V sibiricum var. 
zuccarinii has a maintained moderate level of genetic diversity (mean H E = 0.596, mean N A = 5.6, and 
mean AR = 4.7). Pairwise ^ST indicated significant genetic differentiation between all population pairs 
(0.065-0.309). The semi-natural grasslands in the Aso region have decreased in area over the past de¬ 
cades. Thus, habitat fragmentation may be causing genetic differentiation among populations of V sibiri¬ 
cum var. zuccarinii. Although the populations of this variety have maintained a moderate degree of ge¬ 
netic diversity, the dynamics of gene flow among them should be continuously monitored. 

Key words: conservation genetics, endangered species, grassland, habitat fragmentation, microsatellite, 
vegetative reproduction, Veronicastrum sibiricum var. zuccarinii 


Anthropogenic disturbance is threatening the 
survival of many organisms globally (Pimm et al. 
1995), with many species having gone extinct in 
the past decades (Lawton & May 1995, Regan et 
al. 2001, Smith et al. 1993). Even extant species 
are facing a reduction in population size (Groom- 
bridge 1992). It is estimated that without aggres¬ 
sive conservation measures, about two-thirds of 
all vascular plants may disappear by the end of 
the 21 st century (Pitman & Jorgensen 2002). Ap¬ 
proximately 7,000 vascular plant taxa are recog¬ 
nized in the Japanese flora, one-quarter of which 
is threatened by extinction (Ministry of the Envi¬ 
ronment of Japan 2007, 2012, 2015). Conserva¬ 
tion efforts are therefore urgently needed for 
these threatened plant species. 


Along with habitat destruction and commer¬ 
cial exploitation, stochastic processes, known as 
the “extinction vortex,” have contributed to the 
destabilizing effects of species survival (Gilpin & 
Soule 1986). Some of these stochastic factors are 
genetic, as the decline of genetic variation within 
populations can cause inbreeding depression and 
increase the risk of extinction (Frankham 2005). 
In addition, genetic drift may reduce the genetic 
variation of isolated populations and result in 
strong genetic differentiation among populations 
(Lienert 2004). In small or repeatedly bottle¬ 
necked populations, slightly deleterious muta¬ 
tions are expected to accumulate because of ge¬ 
netic drift, resulting in mutational meltdown and 
a rapid decline toward extinction (Lande 1995, 
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Table 1. The materials and their sources analysed for micro satellite variation 
of Veronicastrum sibiricum var. zuccarinii in Kyushu district, Japan. 


Population 

Locality* 

Altitude 

(m) 

No. of 
samples 

Estimated 
population size 

NA1 

Namino, Aso, Kumamoto 

760 

35 

-100 

NA2 

Namino, Aso, Kumamoto 

780 

35 

-100 

NA3 

Namino, Aso, Kumamoto 

770 

16** 

-100 

NA4 

Namino, Aso, Kumamoto 

770 

33 

-100 

TA1 

Takamori, Kumamoto 

870 

35 

-100 

TA2 

Takamori, Kumamoto 

790 

76 

-500 

KO 

Kokonoe, Oita 

800 

4 

-10 

Total 



234 



*From a view point of conservation, the approximate information of the 
localities is shown. 

**In NA3 population, many ramets grew densely relatively in a small area, 
therefore the collected number of samples were relatively small. 



Fig. 1. Sampling locations of Veronicastrum sibiricum var. zuccarinii in the 
Aso region, Kyushu, Japan. 


Lynch & Gabriel 1990). In threatened plant spe¬ 
cies, the reduction of genetic variation within 
populations may increase both the short-term risk 
of extinction, by reducing reproductive success 
(caused by pollen limitation and inbreeding de¬ 
pression) and the long-term risk of extinction by 
reducing the genetic capacity to respond to envi¬ 
ronmental changes (Frankham et al. 2002, Kery 
et al. 2000). Thus, conservation planning must 
take into account the genetic status of each popu¬ 
lation and individual (Neel & Cummings 2003). 

Veronicastrum sibiricum (L.) Pennell var. 
zuccarinii (Koidz.) H. Hara (Plantaginaceae) is 
an endangered perennial herb of sunny meadows 


in the mountainous areas of southern Korea and 
northern Kyushu, Japan (Yamazaki 1993). Over 
the last few decades, changes in rural land use 
have greatly reduced and fragmented Japan’s 
semi-natural grasslands, and many grassland 
species have become endangered (Tsukada et al. 
2004), V. sibiricum var. zuccarinii being one of 
them. This variety has experienced a rapid popu¬ 
lation decline due to habitat loss resulting from 
abandonment of traditional grassland manage¬ 
ment and conifer plantation establishment, with 
only an estimated seven hundred individuals re¬ 
maining (Environment Agency of Japan 2000). 
Based on the observed rapid decline in habitat 
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remnants and population numbers, var. zuccari¬ 
nii (not the species as a whole) was categorized as 
“endangered” (EN) in the Japanese Red Data 
Book and List (Environment Agency of Japan 
2000, Ministry of the Environment of Japan 
2007). Recently, its status has been revised to 
“vulnerable” (VU; Ministry of the Environment 
of Japan 2012, 2015), because relatively large 
populations have been discovered recently in the 
Aso region. However, currently, population num¬ 
bers are still limited. 

Although Veronicastrum sibiricum var. zuc- 
carinii is at risk of extinction, the population ge¬ 
netics of this variety have not been studied. In the 
present study, we developed and used microsatel¬ 
lite markers to identify the genetic structure and 
diversity of var. zuccarinii, including the fine- 
scale genetic structure within population. Highly 
variable microsatellite markers, which allow for 
accurate estimates of inter- and intra-population- 
al differences, can detect populations with rare 
alleles that may be candidates for targeted con¬ 
servation efforts, as well as populations with low 
diversity that potentially require management as¬ 
sistance (Frankham et al. 2002). We sampled ex¬ 
tant populations from the Aso region. The data 
obtained will help inform the conservation and 
restoration of V sibiricum var. zuccarinii. 

Materials and Methods 


er specimens used in the present study were de¬ 
posited in the herbarium of the Faculty of Sci¬ 
ence, Kumamoto University (KUMA). 

For the TA2 population, we established an 18 
x 16 m quadrat containing 444 individuals (the 
number of ramets). Of these, 66 ramets were ran¬ 
domly selected for collection and the location of 
each was mapped (Fig. 2). The samples were pre¬ 
served as above and used to evaluate the fine- 
scale genetic structure (distribution of genets) of 
the population of Veronicastrum sibiricum var. 
zuccarinii. 

Isolation and characterization of micro satellite 
markers 

The development of microsatellite markers in 
the present study was based on improved tech¬ 
niques for isolating codominant compound mic¬ 
rosatellite markers (Kaneko et al. 2007, Lian et 
al. 2006, Lian & Hogetsu 2002). A restriction 
DNA library was constructed by digesting DNA 
with Sspl or FcoRV, a blunt-end restriction en¬ 
zyme. The restriction fragments were then ligat¬ 
ed using a DNA ligation kit (Takara Bio) with a 
specific blunt-end adaptor (consisting of the 48- 
mer 5’-GTAATACGACTCACTATAGGGCAC- 
GCGTGGTCGACGGCCCGGGCTGGT-3’ along 
with the 8-mer 5’-ACCAGCCC-NH 2 -3’, in which 
the 3’-end was capped with an amino residue). 
Fragments were amplified from the Sspl DNA li- 


Plant sampling 

In July and August 2009-2011, we collected 
mature leaf tissues of Veronicastrum sibiricum 
var. zuccarinii from seven populations (NA1- 
NA4, TA1, TA2, and KO) in the Aso region of 
Kumamoto and Oita prefectures, Kyushu, Japan 
(Table 1, Fig. 1). The collections comprised al¬ 
most all known populations in this region, includ¬ 
ing newly discovered populations (NA3 and 
NA4). Sampling was randomly performed 
throughout each population. A total of 234 indi¬ 
viduals was collected. All samples were dried 
and preserved with silica gel. Population sizes 
were estimated for each population according to 
the number of flowering plants observed. Vouch- 
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Fig. 2. Distribution of 66 ramets within the quadrat estab¬ 
lished in the TA2 population of Veronicastrum sibiri¬ 
cum var. zuccarinii (see, Fig. 1 and Table 1). Same char¬ 
acter indicates plants with identical genotypes; we con¬ 
sidered these plants to be same individuals or clones. 
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Table 2. Characteristics of the five compound micro satellite loci for Veronicastrum sibiricum var. zuccarinii and their variability. 


Locus 

Repeat motif 

Primer sequence (5'-3') 

C(°C) 

Size range (bp) 

A 

H 0 

h e 

Accession No. 

VszOll 

(AC) 6 (AG) 15 

ACACACACACACAGAGAGAGAG 

TTTACCAGTCTTTATCCAACGA 

57 

80-152 

26 

0.604 

0.770 

AB926409 

Vsz016 

(AC) 6 (AG)„ 

ACACACACACACAGAGAGAGAG 

C ATAT C TTT GC AGTAGG AAGC A 

57 

160-190 

7 

0.577 

0.514 

AB926410 

Vsz025 

(AC) 6 (AG) 13 

ACACACACACACAGAGAGAGAG 

AGTTACAGATTTAGCGGCCTAC 

57 

148-208 

24 

0.805 

0.775 

AB926411 

Vsz04 7 

(AC) 6 (AG) 8 

ACACACACACACAGAGAGAGAG 

GTCCTGGCTTTGAGTAACCA 

57 

199-207 

5 

0.479 

0.407 

AB926412 

Vsz226 

(AG) 6 (AC) 10 

AGAGAGAGAGAGACACACACAC 

GAACAATGCTCTGTAACCACAT 

57 

223-227 

7 

0.501 

0.516 

AB926413 

Average 





13.2 

0.593 

0.596 



T a , annealing temperature; A, the number of allele; H 0 , observed heterozygosity; H E , expected heterozygosity 


brary using the compound microsatellite primer 
(AC) 6 (AG) 5 or (TC) 6 (AC) 5 and an adaptor primer 
(5’-CTATAGGGCACGCGTGGT-3’). The ampli¬ 
fied fragments, ranging from 400 bp to 800 bp, 
were then separated on a 1.5% L03 agarose gel 
(Takara Bio) and purified using the QIAquick Gel 
Extraction Kit (Qiagen). The purified DNA frag¬ 
ments were subsequently cloned using the QIA¬ 
GEN PCR Cloning plus Kit (Qiagen) according 
to the manufacturer’s instructions. Briefly, poly¬ 
merase chain reaction (PCR) products were ligat¬ 
ed into the pDrive vector and transformed into 
QIAGEN EZ competent cells. Transformants 
were identified by blue/white screening on LB 
agar plates containing ampicillin, X-gal, and 
IPTG. Cloned fragments containing (AC) 6 (AG) n 
or (TC) 6 (AC) n compound microsatellite sequenc¬ 
es were amplified from the plasmid DNA of the 
positive clones, using M13 forward and reverse 
primers. Based on the sequences obtained, 13 
candidate primers were designed using Primer3 
version 0.4.0 (Rozen & Skaletsky 2000). Six of 
the candidate loci ( VszOll , Vsz016 , Vsz025, 
Vsz047 ’ Vsz205 , and Vsz226) were amplified using 
PCR and showed polymorphic band patterns 
within and among the populations (Table 2). 
Next, we checked the efficiency of the primers 


using 48 plants from four populations (NA1- 
NA4). The linkage disequilibrium test using 
FSTAT version 2.9.3.2 (Goudet 1995) found sig¬ 
nificant genotypic disequilibrium between 
Vsz025 and Vsz205 (P< 0.05); therefore, the 
Vsz205 locus was excluded from the present anal¬ 
ysis. Additionally, significant departure from 
Hardy-Weinberg equilibrium (HWE) after Bon- 
ferroni correction was detected in VszOll , Vsz025 , 
and Vsz226 in the NA1 population, as well as that 
of Vsz047 and Vsz205 in the NA4 population. 
However, no locus showed significant departure 
from HWE in all populations. Thus, a total of five 
loci (excluding Vsz205) was used to assess the ge¬ 
netic diversity of Veronicastrum sibiricum var. 
zuccarinii (Table 2). 

Micro satellite genotyping 

Genomic DNA of all samples was extracted 
from leaves using a modified version of the CTAB 
method (Doyle & Doyle 1987). The five afore¬ 
mentioned microsatellite loci were amplified by 
PCR using a QIAGEN Multiplex PCR Kit (Qia¬ 
gen) (Table 2). Each 5 pL PCR reaction mixture 
contained 0.5 pL template DNA (approx. 1-5 ng), 
1.0 pM primers, and 2.5 pL Master mix. The 
PCR program included initial denaturation at 
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95°C for 15 min followed by 30 cycles of 30 s at 
94°C, 1 min 30 s at 57°C, 1 min at 72°C; these 
were followed by a final extension at 60°C for 30 
min. PCR was performed with the Thermal Cy¬ 
cler 2720 (Applied Biosystems). Sizes of the PCR 
products were estimated using the CEQ™ 8000 
Analysis System with a DNA Size Standard Kit- 
400 (both from Beckman Coulter). 

Data analysis 

We sampled and genotyped 234 plants from 
the seven populations (Table 1). In the following 
analyses, the KO population data were added 
only in the Bayesian clustering analysis (see be¬ 
low) because the sample size was small (four in¬ 
dividuals). To estimate genetic diversity within 
each population, we calculated the mean number 
of alleles per locus (A a ), the expected and ob¬ 
served heterozygosities (H E and H 0 ), the number 
of private alleles ( Pr ) present only in a single pop¬ 
ulation, and the summed number of rare alleles 
with frequencies less than 5% among the total 
population (RA), using GenAlEx version 6.2 
(Peakall & Smouse 2006). Also, allelic richness 
(AR; El Mousadik & Petit 1996) and inbreeding 
coefficients (F IS ) were calculated using FSTAT 
version 2.9.3.2. The latter parameter, F IS , was 
tested with 600 randomizations and by applying 
Bonferroni correction, to determine whether its 
value differed significantly from zero. 

Recent bottlenecks in the populations were 
evaluated using BOTTLENECK ver. 1.2.02 (Piry 
et al. 1999). We simulated equilibrium conditions 
(1,000 replications) assuming the infinite allele 
mutation model (IAM) and the two-phase model 
(TPM) with parameters recommended by Piry et 
al. (1999) (95% single step mutations, variance 
among multiple steps, 12). We used the Wilcoxon 
signed rank test to determine whether a signifi¬ 
cant excess of heterozygosity existed. 

Pairwise Tv (Weir & Cockerham 1984) mea¬ 
sures the degree of genetic differentiation among 
populations and was computed using FSTAT ver¬ 
sion 2.9.3.2. This was tested with 300 permuta¬ 
tions to determine whether it differed significant¬ 
ly from zero. To test whether genetic differentia¬ 
tion among populations followed an isolation-by¬ 


distance (IBD) pattern, a Mantel test (Mantel 
1967) was conducted using GenAlEx version 6.2. 
In this test, the latitude and longitude data were 
used to determine geographic distance among 
populations, and the formula F^fl-F^f) was 
used to determine the genetic distance. Relation¬ 
ships among populations were inferred using a 
neighbor joining (NJ) tree based on the genetic 
distance among populations D A (Nei et al. 1983), 
generated by Population 1.2.30 beta (Langella 
2007). The robustness of each branch was evalu¬ 
ated by bootstrap probability from 1,000 resam¬ 
plings of the genotype dataset. To infer the ge¬ 
netic structure of the populations, Bayesian clus¬ 
tering was implemented using STRUCTURE 
version 2.3.3 (Pritchard et al. 2000). This ap¬ 
proach estimates the number of clusters ( K ) and 
assigns each individual to a cluster. The popula¬ 
tion structure was simulated with values of K = 
1-9 under the LOCPRIOR model with admixture 
model and correlated allele frequencies model (F- 
model). The F-model assumes that all clusters di¬ 
verged from a common ancestral population and 
may have experienced different degrees of genet¬ 
ic drift (Falush et al. 2003). In this model, the 
amount of genetic drift for each cluster is de¬ 
scribed as the F-value. The LOCPRIOR model 
modifies the prior distribution for each individu¬ 
al’s population assignment and allows the propor¬ 
tion of individuals assigned to a particular cluster 
to vary by location (Hubisz et al. 2009). A series 
of 20 independent runs for each value of K rang¬ 
ing from 1 to 9 was performed with U10 6 MCMC 
iterations after U10 5 burn-in periods. To identify 
the appropriate value of K, we compared the 
mean value of the log probability of data (Ln 
P{D)) for each K , according to Pritchard et al. 
( 2000 ). 

Estimation of the fine-scale genetic structure of 
V. sibiricum var. zuccarinii 

To estimate the fine-scale genetic structure 
(distribution of genets) of Veronicastrum sibiri¬ 
cum var. zuccarinii , we genotyped a part of indi¬ 
viduals (66 ramets) collected from the 18x16 m 
quadrat established within the TA2 population 
using the five microsatellite loci and methods de- 
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TABLE 3. Genetic diversity measurements of each population of Veronicastrum sibiricum var. zuccarinii. 


Population 

N 

n g 

n a 

AR 

RA 

Pr 

Ho 

h b 

Fis 

P value of bottleneck analysis 

IAM TPM 

NA1 

35 

16 

5.8 

4.9 

10 

4 

0.613 

0.646 

0.085 

0.031 

0.031 

NA2 

35 

19 

10.2 

7.1 

28 

16 

0.656 

0.701 

0.090 

0.813 

0.813 

NA3 

16 

12 

4.0 

3.7 

5 

1 

0.683 

0.565 

-0.168 

0.156 

1.000 

NA4 

33 

10 

4.2 

4.1 

8 

2 

0.489 

0.563 

0.185 

0.156 

0.625 

TA1 

35 

13 

4.4 

4.1 

8 

1 

0.600 

0.549 

-0.054 

0.156 

0.625 

TA2 

76 

12 

4.8 

4.4 

9 

2 

0.517 

0.556 

0.113 

1.000 

1.000 

KO 

4 

2 

_* 

_* 

_* 

3 


_* 

_* 


_* 

Average** 


13.7 

5.6 

4.7 

10.0 

4.3 

0.593 

0.596 





N, number of samples; Ag, number of genets; Na, mean number of alleles per locus; AR, allelic richness; RA, summed number of rare alleles; 
Pr, summed number of private alleles; Ho, observed heterozygosity; He, expected heterozygosity; Fis, inbreeding coefficient; IAM, infinite 
allele model; TPM, two phase model. 

*In the KO population, the values of genetic diversity within population were not caluculated because the sample size was small. 

**The averages were the value excluding the samples of KO population. 


scribed above. We used GenAlEx version 6.2 to 
calculate the probability of identity, which esti¬ 
mates the average probability that two unrelated 
individuals drawn from the same randomly mat¬ 
ing population will have the same multilocus 
genotype. Because the probability of identity es¬ 
timated by all five loci was very small (3.28x10 ), 
plants with identical multilocus genotypes were 
regarded as clones or genets derived from vegeta¬ 
tive reproduction. We investigated genet distribu¬ 
tion in terms of location within the quadrat and 
individual genotype. 

Observation of underground parts 

To characterize the mode of vegetative repro¬ 
duction, we observed the underground parts of 
Veronicastrum sibiricum var. zuccarinii. In a pri¬ 
vate garden at Takamori-machi, Kumamoto Pref., 
the variety was cultivated in the range of approx¬ 
imately 1 m x 1 m, and the number of the ramets 
with inflorescence was approximately 10. We 
carefully dug up the parts from the stock (approx¬ 
imately 30 cm x 40 cm). After the removal of 
loose soil and rinsing, the underground parts 
were stored in 70% ethanol. 


Results 

Marker loci polymorphisms 

The number of alleles per locus ranged from 4 
for Vsz047 to 26 for VszOll , with an average of 
13.2 (Table 2); a total of 66 alleles were detected. 
Observed heterozygosity ( H 0 ) ranged from 0.479 
for Vsz047 to 0.805 for Vsz025 , with an average of 
0.593. Expected heterozygosity (H E ) ranged from 
0.407 for Vsz047 to 0.775 for Vsz025 , with an av¬ 
erage of 0.596. Nucleotide sequence data report¬ 
ed here are available in the DDBJ/EMBL/Gen- 
Bank databases under the accession numbers 
AB926409-926413. 

Genetic diversity within populations 

We collected 16-76 individuals from each 
population excluding the KO population. Howev¬ 
er, the number of genets ( N G ) that shared geno¬ 
types at all five microsatellite loci ranged from 10 
to 19 (Table 3). The same genets were not detect¬ 
ed among different populations. In each popula¬ 
tion, we excluded clones from subsequent calcu¬ 
lations of genetic diversity (genet-level analyses). 
The mean number of alleles across loci ( N A ) 
ranged from 4.0 for NA3 to 10.2 for NA2, and al¬ 
lelic richness (AR) ranged from 3.7 for NA3 to 7.1 
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Table 4. Pairwise F sx (below diagonal) and geographical distance (km; above diagonal) among 
six populations of Veronicastrum sibiricum var. zuccarinii. 



NAl 

NA2 

NA3 

NA4 

TA1 

TA2 

NAl 


1.0 

3.7 

1.6 

7.1 

12.1 

NA2 

0.089* 


3.7 

2.2 

6.8 

12.3 

NA3 

0.240* 

0.138* 


2.3 

3.6 

8.6 

NA4 

0.171* 

0.137* 

0.218* 


5.9 

10.5 

TA1 

0.065* 

0.105* 

0.284* 

0.275* 


7.0 

TA2 

0.068* 

0.142* 

0.309* 

0.188* 

0.152* 



*Significant differences from zero are indicated (P < 0.05). 


for NA2; mean values were 5.6 and 4.7, respec¬ 
tively (Table 3). Forty-two alleles (64%) were rare 
alleles ( RA ) with frequencies less than 5%. Out of 
a total of 66 alleles, 26 were private alleles ( Pr ), 
which were observed only in a single population. 
The mean number of private alleles per locus 
ranged from 0.2 for NA3 and TA1 to 3.2 for NA2, 
with an average of 0.8. H 0 ranged from 0.489 for 
NA4 to 0.683 for NA3 (mean, 0.593), and H E 
ranged from 0.549 for TA1 to 0.701 for NA2 
(mean, 0.596). F ls ranged from -0.168 for NA3 to 
0.185 for NA4, and no inbreeding coefficient sig¬ 
nificantly differed from zero. In the KO popula¬ 
tion, two genets were detected in four samples 
and three private alleles were detected. 

BOTTLENECK (tested by Wilcoxon’s signed 
rank test) analysis indicated recent population 
bottlenecks in the NA1 population under the I AM 
and TPM (Table 3). There was no excess of het¬ 
erozygosity in any of the other five populations 
under either model. 

Genetic differentiation and structure among pop¬ 
ulations 

Although the geographical distance among 
populations was small (1.0-12.3 km), significant 
genetic differentiation was observed in Veroni¬ 
castrum sibiricum var. zuccarinii. Pairwise F ST 
ranged from 0.065 between the NA1 and TA1 
populations to 0.309 between the NA3 and TA2 
populations (Table 4). F ST values between all pop¬ 
ulation pairs were significantly greater than zero 
(P < 0.05). 


The analysis of genetic structure consistently 
showed that genetic differentiation among popu¬ 
lations did not correspond to geographical distri¬ 
bution. The Mantel test did not indicate a positive 
correlation between geographic distance and ge¬ 
netic distance among populations ( R 2 = 0.005, P = 
0.465) (Fig. 3). The neighbor-joining tree based 
on genetic distances among populations ( D A ) 
showed that all branches had bootstrap probabili¬ 
ties of 50% or less (Fig. 4), indicating that the re¬ 
lationships among the populations were unclear. 

In Bayesian clustering, the mean log probabil¬ 
ity of the data (Ln P(D )) over 20 runs for each 
cluster was highest at K= 7 (Fig. 5). The result of 
a single run at K = 7 is shown in Figure 6. In all 
populations, the cluster with the highest propor¬ 
tion of membership differed (Clusters 1-7). The 
F-value of each cluster ranged from 0 to 0.53, the 
highest value of which occurred in cluster 7. 

Quadrat investigation in the TA2 population 
Analysis of 66 ramets collected from a quadrat 
established in the TA2 population, based on mul¬ 
tilocus genotypes of five microsatellite markers, 
revealed that these ramets consisted of 8 genets 
(Fig. 2). Among the ramets, 56 were members of 
Genet 1, three of Genet 6, and two of Genet 4. The 
remaining ramets were members of Genets 2, 3, 5, 
7, and 8. The longest distance between plants of 
Genet 1 was 1,577 cm (black circles, Fig. 2). 

Observation of underground parts 

The excavated parts included 20-30 ramets of 
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Fig. 3. Relationship between log geographic distance and genetic distance 
between all pairs of populations of Veronicastrum sibiricum var. zucca- 
rinii. 


Fig. 4. Neighbor joining (NJ) tree based on 
genetic distance (DA among popula¬ 
tions of Veronicastrum sibiricum var. 
zuccarinii. Numbers along branches 
are bootstrap values (%) based on 
1,000 replicates. For abbreviations of 
population names, see Table 1. 


Veronicastrum sibiricum var. zuccarinii (Figs. 
7A, B). Rhizomes with horizontal spread were 
observed in the underground parts of these plants 
(Fig. 1C). New shoots for the next year were ob¬ 
served at the base of terrestrial stems (Fig. 7D). 

Discussion 

Moderate genetic diversity in Veronicastrum si¬ 
biricum var. zuccarinii 

To evaluate whether the genetic diversity of 
Veronicastrum sibiricum var. zuccarinii has been 
lost, we compared the mean number of alleles per 
locus (TVa), allelic richness ( AR ), and expected 
heterozygosity (H E ) obtained from microsatellite 
analyses of additional species of endangered 
plants in the semi-natural grasslands of the Aso 
region. For Echinops setifer Iljin (Asteraceae), 
designated as vulnerable (VU) in Japan (Minis¬ 
try of the Environment of Japan 2015), the values 
of N a = 2.1, AR = 1.9, H E = 0.33. Those values 
were determined using eight microsatellite loci 
from 10 populations in the Aso region (Kaneko et 
al. 2009). In Polemonium kiushianum Kitam. 
(Polemoniaceae), also endangered (CR) in Japan, 
the values of N A = 3.3, AR = 2.68, H E = 0.40. 
Those values were determined using 10 microsat¬ 
ellite loci from seven populations in the Aso re¬ 
gion (Yokogawa et al. 2013). In Primula sieboldii 


E. Morren (Primulaceae), designated near threat¬ 
ened (NT) in Japan, the parameters of genetic di¬ 
versity of the Aso population were N A = 9.1 and H E 
= 0.786, based on eight microsatellite loci (Honjo 
et al. 2009). In the vulnerable (VU) species Si- 
lene kiusiana (Makino) H. Ohashi & H. Nakai 
(Caryophyllaceae), five microsatellite loci from 
seven populations in the Aso region revealed the 
following values: N A = 12.0, AR = 10.4, and H E = 
0.791 (Yamasaki et al. 2013). Moreover, Nybom 
(2004), using microsatellite marker data, report¬ 
ed mean H E values in long-lived perennial, nar¬ 
rowly distributed, and outcrossing plants to be 
0.68, 0.56, and 0.65, respectively. In the present 
study of V. sibiricum var. zuccarinii , we deter¬ 
mined the values of N A = 5.6, AR = 4.7, and H E = 
0.596 (Table 3), indicating a moderate degree of 
genetic variation compared to other semi-natural 
grassland species, although the populations of 
var. zuccarinii are currently decreasing in size 
(Environmental Agency of Japan 2000). 

Mode of reproduction is known to influence 
genetic diversity in a species (Hamrick & Godt 
1996, Nybom 2004, Duminil et al. 2007). For ex¬ 
ample, Nybom (2004) suggested that the analyses 
based on microsatellite maker show that long- 
lived and outcrossing taxa retain most of their ge¬ 
netic variability within populations. Veronicas¬ 
trum sibiricum var. zuccarinii is considered to be 
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capable of outcrossing via insect pollination. 
From July through August, V sibiricum var. zuc- 
carinii forms spike-form inflorescences with nu¬ 
merous, crowded flowers (Kitamura et al. 1957, 
Yamazaki 1993). Although a detailed study of the 
breeding system of this species has not been per¬ 
formed, a related species, V japonicum (Nakai) 
T.Yamaz., has hermaphroditic and protandrous 
flowers that are pollinated by bumblebees, flower 
flies, and butterflies (Tanaka 1997, 2009). Similar 
pollinators have also been reported in the related 
North American species, V virginicum (L.) Farw 
(Belt 2011, Levin & Kerster 1969, Vaudo et al. 
2014). Preliminary observations revealed that 
both bumblebees and flower flies visited flowers 
among the populations of V sibiricum var. zucca- 
rinii. Therefore, var. zuccarinii is considered to 
be outbreeding, although it also reproduces veg- 
etatively (for a discussion of vegetative reproduc¬ 
tion, see below). The breeding systems of var. 
zuccarinii may therefore contribute to maintain¬ 
ing moderate genetic diversity within popula¬ 
tions. 

Another possible explanation for the moder¬ 
ate genetic diversity observed in Veronicastrum 
sibiricum var. zuccarinii involves demographic 
processes (population size changes) in the Aso 
region. This semi-natural grassland area has de¬ 
creased from approximately 70,000 ha in 1900 to 
approximately 36,000 ha around 1990 (National 
Parks Association of Japan 1995, Shoji 2006). 
This loss in area has accelerated during the past 
two decades; the grassland area had been reduced 
to 23,000 ha in 2010 (Takahashi 2010). We there¬ 
fore surmise a corresponding decrease in the size 


of populations of V sibiricum var. zuccarinii. In 
small, fragmented populations, random genetic 
drift will reduce genetic diversity (Frankham et 
al. 2002). In fact, the evidence of genetic drift 
was obtained in the Bayesian clustering analysis 
(Fig. 6), and recent bottlenecks were detected in 
the NA1 population (Table 1). However, several 
generations are required before genetic drift has a 
significant impact, especially in long-lived peren¬ 
nial herbaceous plants (Ohara et al. 2006, Tomi- 
matsu & Ohara 2003). In our observations of the 
underground parts of V sibiricum var. zuccarinii , 
we found evidence of vegetative reproduction by 
rhizomes (Figs. 2 & 7), suggesting that the indi¬ 
viduals of V sibiricum may persist for relatively 
long periods. Furthermore, such clonal reproduc¬ 
tion may decrease the rate at which alleles are 
lost, enhance heterozygosity, and thereby in- 



K 

Fig. 5. Distribution of values of log probability of data (Ln 
P(£>)) for each cluster (K) ranging from 1 to 9 calculated 
by Bayesian clustering. 
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Fig. 6. Proportion of membership of clusters in 234 individuals from seven populations of Veronicastrum sibiricum var. zuc¬ 
carinii estimated by Bayesian clustering. 
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Fig. 7. Horizontally spreading rhizomes were observed below ground in Veronicastrum sibiricum var. zuccarinii. A: Habitat 
of observed individuals; B: Upside of stock; C: Closeup of rhizomes; D: New shoots for following year. 


crease the effective population size (Balloux et 
al. 2003, Bengtsson 2003). Hence, despite severe 
habitat reduction, the genetic diversity of V si¬ 
biricum var. zuccarinii may continue to be main¬ 
tained. 

Genetic differentiation and habitat fragmenta¬ 
tion 

The values of pairwise F ST (indicating degrees 
of genetic differentiation among populations) 
were significantly greater than zero between all 
population pairs (mean F ST = 0.17, Table 4). In ad¬ 
dition, Bayesian clustering, according to 
Pritchard et al. (2000), showed that the clusters of 
highest proportion of membership (Clusters 1-7) 
were different in all populations (Fig. 6). Further¬ 
more, both private (Pr) and rare alleles (FA) were 
detected in each population (Table 3). Therefore, 
we consider genetic differentiation to have oc¬ 
curred among populations of Veronicastrum si¬ 
biricum var. zuccarinii. 

Over the past several decades, as mentioned 


above, the habitat of Veronicastrum sibiricum 
var. zuccarinii (and other grassland plants) has 
decreased due to the decline in semi-natural 
grassland areas in the Aso region (National Parks 
Association of Japan 1995). In plants that cannot 
achieve long-distance gene flow, it is known that 
greater distances between populations because of 
habitat fragmentation decrease the likelihood of 
gene flow (Primack & Miao 1992, van Dorp et al. 
1997, Cain et al. 2000, Soons et al. 2005). Fur¬ 
thermore, insect-pollinated species with gravity- 
based seed dispersal are believed to have de¬ 
creased gene flow between populations when 
compared with wind-pollinated plants and those 
with animal- or wind-dispersed seeds (Loveless 
& Hamrick 1984). As mentioned above, the spe¬ 
cies of Veronicastrum are presumed to be polli¬ 
nated by insects (bumblebees, flower flies, and 
butterflies) (Tanaka 1997, 2009, Belt 2011, Vaudo 
et al. 2014). Furthermore, V. sibiricum forms cap¬ 
sules with many small (approximately 0.6 mm 
long), hemispherical seeds (Yamazaki 1993). We 
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surmise that the seeds are mainly dispersed by 
gravity because the seeds and capsules lack 
wind-dispersal appendages. Veronicastrum si- 
biricum var. zuccarinii is most likely pollinated 
by insects and produces gravity-dispersed seeds, 
making it more susceptible to the effects of habi¬ 
tat fragmentation. 

Although we inferred genetic differentiation 
among populations of Veronicastrum sibiricum 
var. zuccarinii , the Mantel test did not reveal a 
pattern of isolation by distance (IBD) (Fig. 3). In 
the case of inter-population differentiation with¬ 
out IBD, it is estimated that the effects of genetic 
drift in each population are greater than those of 
gene flow (Hutchison and Templeton 1999). In 
fact, the Bayesian structure analysis (Fig. 6), in¬ 
dicating the degree of genetic drift for each clus¬ 
ter (Falush et al. 2003), showed relatively high up¬ 
values, except for Cluster 2 (0.18-0.53), and re¬ 
cent bottlenecks were also detected in the NA1 
population (Table 1). Therefore, it is considered 
that populations of V sibiricum var. zuccarinii 
are greatly affected by genetic drift because of 
habitat fragmentation and differences in genetic 
composition among the population are increas¬ 
ing. 

Conservation perspective 

We investigated the genetic diversity of Ve¬ 
ronicastrum sibiricum var. zuccarinii from seven 
populations, which constitute nearly all the 
known populations of this variety in the Aso re¬ 
gion of Japan. Our study revealed that each popu¬ 
lation has maintained moderate genetic diversity. 
However, we found evidence of significant genet¬ 
ic differentiation among populations (Table 4). In 
addition, populations of V sibiricum var. zuccari¬ 
nii in the Aso region are evidently small (estimat¬ 
ed roughly as 100-500 ramets/population) (Table 
1). Small, fragmented populations are subject to 
extinction or substantial reductions in population 
size due to stochastic events (Lande 1993, Mel¬ 
bourne & Hastings 2008), which may lead to ge¬ 
netic erosion depending on the degree to which 
the effective population size is reduced. Further¬ 
more, in most populations, the number of genets 
is smaller than the number of ramets (Table 3) 


and we found individuals in the TA2 population 
to be actively capable of vegetative reproduction 
by rhizome (Figs. 2 & 7). By enhancing the vari¬ 
ance in genet size and reducing the opportunities 
for recombination, such clonal reproduction may 
also reduce the effective population size and con¬ 
sequently enhance the magnitude of genetic drift 
(Orive 1993, Campbell & Husband 2005). There¬ 
fore, we recommend that continuous monitoring 
be carried out on the population dynamics of V 
sibiricum var. zuccarinii. 

In the present analysis, the NA2 and NA1 
populations showed higher genetic diversity (N A , 
AR , RA , Pr, and H\) than the other populations 
(Table 3). Furthermore, the degree of genetic drift 
in the population was estimated to be relatively 
small (F-values = 0 and 0.18) in the Bayesian 
structure analysis (Fig. 6). Although we do not 
know of any specific factors with respect to the 
genetic status of the populations, both popula¬ 
tions have been protected as a semi-natural grass¬ 
land for many years. The grasslands have been 
managed by burning and mowing by the Aso 
Green Stock Foundation, an organization devoted 
to grassland regeneration (Yamauchi & Taka- 
hashi 2002). Accordingly, in these populations, 
suitable genetic traits may be maintained longer 
than in other populations. From the viewpoint of 
conservation genetics, the NA2 and NA1 popula¬ 
tions are very important and should be preferen¬ 
tially conserved, although we should conserve as 
many as possible of the populations, because our 
analyses showed V sibiricum var. zuccarinii ex¬ 
hibits significant genetic differentiation among 
populations (Tables 3 & 4, Fig. 6). 

We estimated the genetic diversity of Veroni¬ 
castrum sibiricum var. zuccarinii based on the 
five microsatellite markers used in the present 
study. In general, 10 or more markers are fre¬ 
quently used to evaluate the genetic diversity and 
structure of a population. It is known that esti¬ 
mates of average heterozygosity and genetic dis¬ 
tance are affected by the number of loci analyzed 
(Gorman & Renzi 1979, Nei 1978, Nei & Roy- 
choudhury 1974, Nybom 2004). Therefore, with 
respect to these parameters, the results should be 
re-evaluated by increasing the number of markers 
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that are analyzed. 
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